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The crossover from Edwards- Wilkinson (s = 0) to KPZ (s > 0) type growth is 
studied for the BCSOS model. We calculate the exact numerical values for the k = 
and 271 /N massgap for < 18 using the master equation. We predict the structure 
of the crossover scaling function and confirm numerically that mo — 4(7r/A^)^[l + 
3w2(s)Ar/(27r2)]0-5 and mi ~ 2(7r/A^)2[l + u'^{s)N/T:'^f-^ , with u{l) = 1.03596967. 
KPZ type growth is equivalent to a phase transition in meso-scopic metallic rings 
where attractive interactions destroy the persistent current; and to endpoints of 
facet-ridges in equilibrium crystal shapes. 
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A large amount of theoretical effort has been devoted in recent years to establish 
and classify the scaling properties of dynamic processes such as those at growing 
interfaces. Several dynamic universality classes have emerged. One of them is 
the so-called KPZ universality class, named after the non-linear Langevin equation 
studied by Kardar, Parisi, and Zhang [1]. Older examples are Edwards- Wilkinson 
(EW) growth [2] and directed percolation [3] . They are distinguished by the value 
of the dynamic critical exponent z; characteristic lengths and times scale as i ~ P. 
Most of the evidence is numerical in nature, in particular from Monte Carlo (MC) 
studies [2], supplemented by a few exactly soluble models [4] and field theoretical 
renormalization studies [1,2]. 

Our understanding of dynamic scale invariance is still nowhere near that of 
equilibrium critical phenomena. In particular the comparison with two dimensional 
(2D) critical phenomena is relevant. ID dynamic processes and 2D equilibrium sta- 
tistical mechanics are very much alike in the master equation and transfer matrix 
formulation. Time becomes the second spatial coordinate. Most 2D equilibrium 
critical phenomena obey conformal invariance [5]. They all belong to the same 
conformal "dynamic" universality class (2; = 1; rotational invariance). The ground 
state properties of ID quantum systems are related as well. Luttingcr liquids (rela- 
tivistic fermions) are also part of the conformal z = 1 dynamic universality class [6] . 
Non-relativistic fermions have a dynamic exponent z = 2. Anisotropic scaling, with 
2; 1, occurs in 2D critical phenomena at e.g. Pokrovksy-Talapov (PT) transitions 
{z = 2) and the elusive Lifshitz points [6] . One prospect is that concepts similar to 
conformal invariance might apply to dynamic universality classes in general. 

In this letter we present a detailed numerical study of the finite size scaling 
(FSS) crossover function from EW {z = 2) to KPZ {z = 1.5) type growth in the 
body-centered solid-on-solid (BCSOS) growth model. It is important to know the 
details of the crossover scaling between dynamic universality classes. The following 
simple example illustrates this. 

Consider the ID free fermion model. The fermions are non-relativistic when 
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the chemical potential coincides with the bottom of the energy band E{k) ~ k'^. 
The momenta are spaced as Ak = 2tt/N. So the massgap scales as m ~ (27r/A^)^ 
and zq = 2 because It ~ . The fermions are relativistic for /c^ > 0. The 
massgap scales as m ~ i^Tikp/N ^ and zi = 1. The crossover scaling function m = 
{kp + 2tt/NY - kjr ^ A-Kkp/N + {27t/N)'^ has two important properties. The 
crossover exponent j/c determines the scaling of dm/dkp ~ JSf-^o+Vc = q. 

?/c is equal to the change in the dynamic exponent, Uc = zq — zi = 1. The FSS 
amplitude of the massgap is not a constant in the metal phase, but a universal 
number multiplied with the rapidity 2k p- The metal phase belongs to the conformal 
z = 1 dynamic universality class. The rapidity represents the lattice anisotropy in 
2D critical phenomena language. The crossover operator is the stress tensor which 
lies at the heart of conformal field theory. We expect that crossover scaling functions 
between dynamic universality classes have the same generic assymptotic form 

m ~ A/N^° + Bu{s) /N'' , (1) 

with A and B universal constants and u{s) the scaling field associated with the 
crossover operator. 

Surface roughness is characterized by the scaling of the height-height correlation 
function, G{r,t) =< [h{ro + r) -h{ro)]'^ >~ lP'°'G{h~^r,¥t~^). In ID the stationary 
state has typically a finite correlation length beyond which the steps are disordered, 
and the probability to go up or down along the interface becomes random. At large 
length scales G{r) behaves like the root mean-square displacement of a random walk, 
G{r) ~ r. So typically a = 0.5. Stationary states with interesting correlations like 
those in 2D driven diffusive systems [7] have not been found (yet) in ID growth 
dynamics, z characterizes the rate at which the stationary state is approached. A 
disturbance with characteristic size ^ spreads-out as ^ ~ t^. 

In the BCSOS model nearest neighbour columns must differ in height by one 
unit, dh = ±1. This describes a body-centered type stacking of molecules, or a 
zig-zag type stacking of bricks. We represent the surface configuration in terms of 
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the steps, = ±1. The dynamic rule is as follows. Select one of the columns at 
random; if = —1 and 5'^_|_i = 1 (a local valley) a brick is adsorbed with probabil- 
ity p; if S"^ = 1 and S^_^_i = —1 (a local hill top) evaporation occurs with probability 
q. P{{S^}) is the probability of micro state {S^}. |^')^ = J2 P{{Sn})\{Sn}) is the 
state vector. The time evolution operator T of the discrete-time master equation 
|*)t+a = ^|*)r has the following structure 

T = 1- -eN ^ ^[1 - XS^S^_^_i - -(5'+<S'^^.i + S^+i) 

n 

~2^i^n^n+l ~ ^n+l) + f^^ril (2) 

with A = l/AT, s = {p — q)/{p + q), and periodic boundary conditions, <S'^+jv = S^. 
Without loss of generality time can be rescaled to e = p + q = 1. The parameters 
A = 1 and /i = are introduced for later convenience. Stochastic processes preserve 
probability. So the state where all microstates have the same probability, \D) = 
S K'^n})) is the left eigenvector for the largest eigenvalue, Aq = 1. 

The BCSOS model has more symmetry than generic KPZ growth models: 
T^iPiO) = ^(9)P)^- The left eigenvectors (L| of T{p,q) are identical to the right 
eigenvectors of T{q^p). T{q,p) and T{p,q) are equivalent by particle-hole symme- 
try. Therefore \R) = [Y[Sn] \L), and \D) is both the left and the right eigenvector 
for Aq. The stationary state is completely disordered; a = 1/2 for all s. The BC- 
SOS model is used widely to study KPZ type growth [2]. This special symmetry 
casts some doubts on its role as a generic KPZ growth model. This becomes more 
apparent in the crossover scaling behaviour. 

It is important to point out how this model fits into the general picture sketched 
above. T has the generic form, T = exp(— 7Y). Ti. is non- Hermitian, except at s = 
where the model reduces to the spin- 1/2 XXZ chain [6]. The EW point at s = 
maps exactly onto the ferromagnetic Heisenberg point, A = 1. In the KPZ equation 
EW growth corresponds to the point where the non-linear term vanishes and the 
dynamics reduces to a stochastic diffusion equation. Indeed, the exact solution of 
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the XXZ chain at A = 1 is diffusive. The equation of motion for the spin-spin 
correlation closes and is a deterministic diflFusion equation [8]. We checked that 
also the equations of motion for the n-point spin correlation functions close within 
themselves and are diffusion type equations. 

The XXZ model is equivalent to a ID spinless interacting fermion problem [6] . 
Growth, s, introduces a preferred hopping direction. It creates a persistent current 
around the ring, similar to those in metallic rings in perpendicular magnetic fields. 
Recent research is focused on coherence effects in the presence of elastic scattering 
from random impurities [9]. Eq.(2) does not include impurities. It describes a 
phase transition due to interactions. Weak attraction between the fermions, small 
A, only renormalizes the scaling indices due to the marginal operator in the chiral 
Luttinger liquid. At A > 1 the fermions coalesce into a macroscopic bound state (the 
ferroelectric ground state of the XXZ chain). This destroys the persistent current, 
because in eq.(2) only single fermions can hop. (In metallic rings the macroscopic 
droplets of fermions might still be able to move). The transition takes place at 
the KPZ line, at A = 1. At this transition the fermion excitations do not behave 
metallic, z = 1, but are described by the unconventional dynamic exponent z = 1.5. 
It will be interesting to investigate how elastic impurities modify this. 

The spin-1/2 XXZ chain is related to the 6- vertex model for 2D equilibrium 
(anti-) ferroelectrics [10]. The equivalence becomes exact when in the dynamic rule 
the coUumns are chosen sequentually instead of randomly, first all even and then all 
odd coUumns [8]. The modified rule reduces to the old one in the time continuum 
limit, p ^ and q ^ with s constant, /i and s are the two components of the 
electric field. 

The 6-vertex model is equivalent to the BCSOS model for equilibrium crystal 
surfaces. The anti-ferroelectric side of its phase diagram describes surface rough- 
ening transitions [10]. The ferroelectric side is less well known. It describes the 
crystal shape close to a facet ridge [11,12]. A plays the role of inverse temperature. 
At A > 1 two facets meet at a facet ridge; see Fig. 1. The growth parameter s tilts 
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the surface in the time-hke direction, the direction along the facet ridge. The chem- 
ical potential /x tilts the surface in the spatial direction, the direction _L to the facet 
ridge. The interactions in this direction are anti-ferrmagnetic and favour facetting. 
The KPZ points are the endpoints of the facet ridge where the two facets become 
separated by two PT transition lines with rough rounded surface in-between; see 
Fig. 1. The KPZ line represents the change in location with temperature of these 
facet ridge endpoints. The facet ridge shortens with temperature until it vanishes 
at the EW critical point. At PT transition the surface rounds smoothly. At facet 
ridge endpoints it changes discontinuously, because at KPZ points the surface has 
a finite slope in the time-like direction equal to the growth rate Vg = s/4. The 
dynamic exponent z = 1.5 implies unusual anisotropic scaling behaviour. There is 
one caveat to this identification of KPZ growth with facet ridge endpoints. More 
general growth models such as the RSOS model lack the self-adjointness property 
of eq.(2). Moreover, at the facet ridge endpoints in the BCSOS model the transfer 
matrix happens to be stochastic; we see no reason for this to be true in general. It 
will be interesting to investigate how these aspects affect the scaling properties. 

The eigenvalues of the transfer matrix can be classified according to wave 
numbers k, reflecting translational invariance. The largest eigenvalue, Aq = 1, 
is located in the k = sector. Define two massgaps: mo = Ao(0) — Ai(0) and 
^1 — ^o(O) — Ao(l). Ai(0) is the next largest eigenvalue in the k = sector, and 
Ao(l) the largest eigenvalue in the k = ^tt/N sector. We calculated the exact nu- 
merical values of mo and mi for even system sizes AT < 18. The dynamic exponent 
follows from the FSS of the massgaps, m ~ N~^. The spatial exponent follows 
independently from the interface width, W ~ A^", in the stationary state (the right 
eigenvector of Ao(0)). 

FSS approximants for z are shown in Fig. 2 at various values of s. z{N, N+2) = 
ln{mi{N)/mi{N + 2))/\n{{N + 2)/N). These numbers are accurate to better than 
12 decimal places, unlike results from e.g. MC simulations. Therefore it is possible 
to take into account the leading corrections to scaling. We do this by constructing 
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a cascade of 1/N extrapolations for the z{N,N + 2) at successive values of N. At 
s = 1 the convergence is very stable, and leads to 2; = 1.50 ± 0.05. This is the same 
level of accuracy as obtained by MC simulations, but with substantially less effort 
and much smaller system sizes. 

The entire line < s < 1 belongs to the KPZ universality class. The dynamic 
exponent must be equal to 2; = 1.5 for all s 7^ 0, but the convergence is only good 
for s ~ 1. For s ^ 0, the z(N, N + 2) move towards the EW value 2; = 2 at small N 
before they bend-down towards the KPZ value z = 1.5. This reflects the crossover 
scaling behaviour from the nearby EW point. 

The crossover scaling function eq.(l) does not apply to the BCSOS model, due 
to the special symmetry discussed above, which implies that the massgaps are even 
functions of s. Their first derivative vanishes. Eq.(l) is likely to apply to the more 
generic growth models, but this needs still to be tested. We propose the modified 
form for the BCSOS model: 

m = ^[l + Bu'{s)Nr-\ (3) 

Also in this expression j/s is equal to the change in the dynamic exponent, and the 
amplitude at s is equal to a universal number multiplied with u{s). 

The crossover scaling exponent follows from the second derivatives of the mass- 
gaps with respect to s, m" ~ N-^+'^Vs ^ We find ya = 0.50 ± 0.05; consistent with 
Us = 0.5 from power counting at the EW point in the KPZ equation. 

At s = the k = 27r/N mass gap is exactly equal to mi (0) = 2 sin^ (t^/N) for all 
N. Eq.(2) reduces to the ferromagnetic Heisenberg model. Spin rotations leave T 
invariant, but mix-up the sectors of different magnetization (slopes of the surface). 
Since Xo{k) is non-degenerate in each sector it follows that nii is independent of 
the tilt of the surface. The equation for mi (0) is trivially true in one particle sector 
with surface tilt 1 — 1 /N. 

The k = massgap at s = is more complex. We find numerically: mo(0) = 
4sin2(7r/iV)[l + Jsin(7r/A^) + (0.306 ± 0.003) sin2(7r/A^) + ...]. The exact result for 
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mi(0) suggests that sin(7r/A'") is the natural FSS parameter. Only FSS corrections 
with integer powers of sin(7r/A^) are to be expected, because all n-point correlation 
functions at s = obey diffusion type equations. We find no numerical evidence 
otherwise. We expect amplitudes to be simple combinations of n and integers. The 
leading correction to scaling amplitude of mo(0), 0.635 ± 0.005, is close enough to 
2/7r to guess this to be the exact value. Similarly, the second derivatives of the 
massgaps at s = scale numerically as: m'/(0) = |- sin(7r/A^) — sin^(7r/A^) + (0.25 ± 
0.02)sin^(7r/A^) + ...] and m[)'(0) = f sin(7r/A^) - (1.06±0.01) sin2(7r/A^) + ...]. These 
results yield the amplitudes in eq.(3): 

mi{s,N) = 2sin2(7r/iV)[l + -w(s)7sin(7r/iV)]°-^ 

TT 

mo{s,N) = 4sm'^{n/N)[l + —u{s)^ / sm{n /N)]^-^ (4) 

27r 

At s = 1 the leading amplitude of the k = 2tt /N massgap is known exactly from 
the Bethe ansatz solution, mi ~ 1.1689666(7r/A^)i [4]. We use this in our scaling 
analysis at s = 1: mi(l) = 1.1689666 sini (7r/iV)[l + (0.531 ± 0.002) sin(7r/7V) + 
...]. This gives information about the scaling field u{s). The renormalization is 
small. The amplitude predicted by eq.(4) assuming tt = s is only 7% too small. 
Our knowledge of the non- linear scaling field can be summarized as i((s) = s(l + 
0.0359697s^/(s^)), with /(s^) an unknown function, except for /(I) = 1. 

The amplitude of the /c = massgap at s=l is not known from the Bethe Ansatz 
solution [4]. The two massgaps must depend on u{s) in the same way. Therefore 
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eq.(4) predicts the exact value, mo(l) 2.8633717 sin2 (tt/A/"). This agrees well 
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with a free numerical estimate, mo(l) ^ (2.87 ± 0.01) sins (tt/A^). Moreover, the 
leading corrections to scaling amplitude is very stable when we assume this value 
to be exact: mo(l) ~ 2.8633717 sini (7r/7V)[l + (0.940 ± 0.005) sin(7r/7V) + ...]. 

It is tempting to visualize the low lying excitations in terms of quasi-particles 
with a dispersion relation E{k) ~ . The amplitude of the massgap mi ~ AN~^ 
yields the slope of such a dispersion relation. In this picture, mo must have an 
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amplitude twice a big as mi, because the lowest k = excitations involves the 
creation of two quasi particles of opposite momentum. This is true at the EW 
point, but not for KPZ type growth where the amplitudes differ by more than a 
factor 2. So the large length scale behaviour of KPZ growth can not be described 
by a free field theory this simple in structure. 

Fig. 3 shows the convergence of the ratio between the numerical values of mi (s) 
and those predicted by eq.(4), as function of s. The leading corrections to scaling 
are incorporated. The representation of corrections to scaling is not unique because 
their origins cannot be distinguished numerically. We choose to represent them as: 

mi = 2sin2(^)[l + - h{s)s' + ...]°-^ 

^N^^ 7rsm(7r/iV) ^ ^ ^ 

with bi{s) — ^ + 0.137s^ from the corrections to scaling in m'/(0) and mi(l). The 
ratio is a sensitive test because both quantities vanish in the thermodynamic limit. 
The k = 27r/N massgap ratio converges to 1 within 0.1% for all s. Similarly, the 
k = massgap ratio converges to within 0.5%. The quality of the convergence 
varies slightly with s. This reflects the uncertainty in the scaling field u{s). The 
shape of the curves is sensitive to the precise choice for /(s^). /(s^) = 1 can be 
ruled out, because our numerical estimates of the amplitudes at s = 0.6 is slightly 
but significantly smaller than predicted by f{s^) = 1. We investigated the simple 
form /(s^) = (1 — a)s^ — a for various values of a. /(s^) = s^, shown in Fig. 3, gives 
the smoothest curves, but the differences are marginal. The convergence of both 
ratios stays within the limits quoted above for any of these choices. 

We conclude that the asymptotic crossover scaling function eq.(3) is correct. 
The agreement with the exact numerical results is within 7% over the entire range 
< s < 1 when s is assumed to be a pure scaling field. The agreement improves to 
within 0.5% when the renormalization of the growth parameter is incorporated. 

This research is supported by NSF grant DMR-9205125. 
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Figure Captions 

Fig.l: Schematic temperature evolution of a facet ridge in an equilibrium crystal. 

KPZ growth coincides with the facet edge endpoint (a) and EW growth with 

the point where the facet ridge vanishes (b). 
Fig. 2: Finite size scaling approximants of the dynamic exponent z at various values 

of s from the k = ^tt/N massgap. 

Fig. 3: The ratio between the numerical values of the k = 2n/N massgap and those of 
the proposed crossover scaling function eq.(4) for system sizes A'" = 10 — 18. 
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